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ABSTRACT 



In Loeb & Rybicki (1999; paper I) it was shown that before reionization, the scattering of 
Lya photons from a cosmological source might lead to a fairly compact (~ 15") Lya halo around 
the source. Observations of such halos could constrain the properties of the neutral intergalactic 
medium (IGM) , and in particular yield the cosmological density parameters of baryons and matter 
on scales where the Hubble flow is unperturbed. Paper I did not treat the polarization of this 
scattered radiation, but did suggest that the degree of such polarization might be large. In this 
Letter we report on improved calculations for these Lya halos, now accounting for the polarization 
of the radiation field. The polarization is linear and is oriented tangentially to the projected 
displacement from the center of the source. The degree of polarization is found to be 14% at 
the core radius, where the intensity has fallen to half of the central value. It rises to 32% and 
45% at the radii where the intensity has fallen to one-tenth and one-hundredth of the central 
intensity, respectively. At larger radii the degree of polarization rises further, asymptotically to 
60%. Such high values of polarization should be easily observable and provide a clear signature of 
the phenomenon of Lya halos surrounding sources prior to reionization. 

Subject headings: cosmology: theory - line: profiles 



1. Introduction 



High-redshift galaxies are detected at present out to z <J 5.6, and are found to be strong 
Lya emitters (Dey et al. 1998; Hu, Cowie, & McMahon 1998; Spinrad et al. 1998; Weymann et 
al. 1998). Popular cosmological models predict that at somewhat higher redshifts, z ~ 10, the 
hydrogen in the intergalactic medium (IGM) was transformed from being predominantly neutral 
to being ionized due to the UV radiation emitted by the first stars and mini-quasars (see, e.g. 
Gnedin &: Ostriker 1997; Haiman & Loeb 1998a, b). Prior to this epoch of reionization, the neutral 
IGM was highly opaque to resonant Lya photons. Hence, the Lya photons emitted by early 
galaxies were scattered in their vicinity by the surrounding IGM. In a previous paper (Loeb & 
Rybicki 1999, paper I), we have shown that this intergalactic scattering results in compact (~ 15") 
halos of Lya light around such sources. The scattered photons compose a line of a universal 
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shape, which is broadened and redshifted by ~ 10 3 km s^ 1 relative to the source. The detection of 
these intergalactic Lya halos could provide a unique tool for probing the neutral IGM before and 
during the epoch of reionization. In addition, we have found that observations of the Lya intensity 
profile on scales where the Hubble flow is only weakly perturbed, could constrain the cosmological 
density parameters of baryons (Ob) and matter (Qm)- 

Paper I had suggested, but not demonstrated, that the scattered Lya light would be highly 
polarized. The polarization signal is important in that it provides an unambiguous signature of 
the scattering nature of the diffuse Lya halos around high-redshift galaxies. Due to the spherical 
symmetry of the scattering geometry, the linear polarization of the scattered radiation is expected 
to be oriented tangentially relative to the projected displacement from the center of the source. 

In this Letter we report on a detailed calculation of the polarization properties of scattered 
Lya halos. In §2 we describe the Monte-Carlo approach employed in this calculation, and in §3 
we describe our numerical results. Finally, §4 summarizes the implications of these results. 

2. Polarized Monte Carlo Method 

Let us first discuss the atomic scattering process for the Lya line. We note that the hydrogen 
Lya line at vq = 2.466 x 10 15 Hz is actually a doublet consisting of the two fine-structure lines, 
2 5i/2 ~ 2 P\/2 an< ^ 2 ^i/2 ~ 2 -^3/2> separated by 1.1 x 10 10 Hz. Fortunately, as shown in paper I, 
the regime of interest to us involves frequency shifts from these line centers of order z/* s=s 10 13 
Hz, which are much larger than the separation of the lines. In this regime, quantum-mechanical 
interference between the two lines acts in such a way as to give a scattering behavior identical to 
that of a classical oscillator, that is, the same as pure Rayleigh scattering (Stenflo 1980). (Using an 
incoherent superposition of the results of Hamilton [1947] for the two lines, one would incorrectly 
conclude that only one-third of the scattering is polarized.) 

We shall now describe the modifications of the Monte Carlo method of paper I necessary to 
treat polarization. Monte Carlo methods for polarized radiative transfer are often formulated 
using "photons" that are actually groups of photons with specified Stokes parameters (see, e.g., 
Whitney 1991; Code & Whitney 1995). However, for the present case we found that a more 
convenient description of polarization was to use individual photons, each with a definite state 
of 100% linear polarization, a description previously used by Angel (1969). (There is no need to 
consider circular polarization here, since the central source is assumed to be unpolarized, and 
Rayleigh scattering cannot generate circular polarization, except from circular polarization.) If 
the direction of the photon is given by the unit vector n, then its polarization is defined by a real 
unit vector e, with n • e = 0. In this formulation, the observed Stokes parameters result from the 
statistics of binning together multiple, independent photons. 

With this description of photons, the polarized Monte Carlo method involves much the same 
steps as the unpolarized version of paper I. One difference is in the handling of the polarized 
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Rayleigh scattering process, which is done as follows: The angular distribution of the scattered 
photon has a probability distribution per solid angle proportional to sin 2 0, where is the angle 
between the scattered photon and the polarization vector of the incident photon. This is simulated 
using a rejection technique: We choose a random unit vector n' (uniform in solid angle) and a 
uniform random deviate R on (0, 1), and test whether R < 1 — (e ■ n') 2 ; if not, we start again 
with new random choices for n' and R; the process is repeated until the test is passed, and then 
n' is taken as the new photon direction. The new polarization vector e' is determined by finding 
the normalized projection of the old polarization vector e onto the plane normal to n', that is, 
e' = g/|g|, where g = e — (e • n')n' . As in paper I, an individual photon is followed through a 
number of scattering events until it escapes. 

The only remaining question is how to characterize the polarization of the escaped, observed 
radiation in the plane of the sky as a function of impact parameter p. From symmetry, we know 
that this radiation can be characterized by the intensities parallel to the projected radius vector, 
Ii, and perpendicular to it, I r . If \ is the angle between the photon's polarization vector and the 
projected radius vector, then it contributes to the appropriate histogram bins a fractional photon 
number cos 2 \ to and sin 2 \ to I r (these are the squares of the components of the polarization 
vector). With appropriate normalizations (see paper I), these histograms determine the two 
intensities and I r . In terms of these the degree of polarization is II = \Ii — I r \/(Ii + I r ). 

An alternative way of stating the results is in terms of the Stokes parameters / and Q, which 
are related to the above intensities by / = + I r and Q = 1\ — I r . The degree of polarization is 
II = \Q\/I. These Stokes parameters can also be found directly in the Monte Carlo method by 
binning with fractional photon numbers 1 for /, and cos 2 x — sin 2 \ = cos 2% for Q. 

The above polarized Monte Carlo method was tested by solving the classical Milne 
problem for a Rayleigh scattering atmosphere (Chandrasekhar 1950; §68). The results for the 
emergent intensities and I r agreed, to within statistical errors, with precise analytical results 
(Chandrasekhar 1950; table XXIV, p. 248). 



3. Results 



As in paper I, we use rescaled variables (denoted by tildes), allowing one single solution to 
apply to all physical cases. In particular, we normalize frequencies by, 
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and distances by, 
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Hz, (1) 
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where Qb, and Q\ are the density parameters of baryons, matter, and vacuum, respectively; z s 
is the source redshift; and ho is the Hubble constant in units of 100 km s _1 Mpc -1 . The radiation 
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intensity is normalized by J* = 7V a /(r*z/*), where N a is the steady emission rate of Lya photons 
by the source. 

The Monte Carlo method was used to follow the scattering of 10 8 photons, a sufficiently large 
number to provide reasonable statistical accuracy (except at very small impact parameters). The 
profile of the frequency-integrated total Stokes intensity / versus impact parameter p is given in 
Figure 1 (upper solid curve). The corresponding result from paper I, in which the scattering was 
approximated as unpolarized and isotropic, is given as the dashed curve. The new curve is seen 
to be slightly more centrally concentrated, with a central intensity about 25% higher than that of 
paper I. 

The "core" radius, where the total Stokes intensity / has fallen to half its central value, is 
at p = 0.070. The radii where the intensity falls to one-tenth and one-hundredth of the central 
intensity are p = 0.25 and 0.80, respectively. At large impact parameters, / falls approximately as 
p~ 3 . 

The polarized intensities Ii and I r are also plotted in Figure 1. At the center of the observed 
disk, p = 0, both of these intensities are equal to 1/2 by symmetry. For all other values of impact 
parameter, we note that I r exceeds Ii everywhere, and asymptotically by a factor of four. Thus the 
radiation is strongly linearly polarized with orientation tangential to the projected radius vector. 

The lower panel of Figure 1 shows the degree of polarization II versus impact parameter. 
This parameter rises monotonically from zero at p = (as required by symmetry), to about 14% 
at the core radius, and to 32% and 45% at the one-tenth and one-hundredth intensity points, 
respectively. At still larger impact parameters, the polarization rises even further to an asymptotic 
value of 60%. However, since the intensities are falling so rapidly, polarizations of that asymptotic 
magnitude may not be detectable in practice. 

Some heuristic insight into why these polarizations are so high can be found by comparison 
to the problem of Rayleigh scattering of radiation from a point source surrounded by an optically 
thin scattering medium with power-law density profile N(r) oc r~ n , first treated by Schuster (1879) 
in the context of Thomson scattering in the solar corona. Schuster showed that the total intensity 
is a power law in impact parameter, / oc p~( n+1 \ and the degree of polarization is (n + l)/(n + 3). 
To relate the asymptotic results of our problem to those of the Schuster problem, we first note 
that the radiation at large impact parameters is dominated by the scattering of photons that have 
taken one very large step from much nearer the source, so they are travelling almost radially, as 
they would from a point source. The second observation is that the Schuster result really applies 
more generally to a power law dependence of the product of the scattering cross section a(r) times 
the density of the form a(r)N(r) oc r~ n . This is because only the optical depth along the radial 
direction is relevant to Schuster's derivation. In the situation treated by Schuster, the scattering 
cross section was constant. In our case the density is constant, but the scattering cross section 
seen by these radial photons decreases inversely as the square of the distance, since the frequency 
displacement is redshifting linearly with distance, and the line profile varies inversely as the 
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square of frequency. Therefore we should compare with the Schuster results for n = 2. This gives 
an inverse third power behavior for the intensity and 3/5 = 60% for the degree of polarization, 
exactly as seen in our asymptotic results. 

The image of the halo taken in unpolarized light will have circular symmetry on the sky. A 
contour plot of such an image is shown in the left panel of Figure 2. The ten contours are spaced 
by one-half magnitude (0.2 dex), so that the outer contour represents an intensity 10 -2 times the 
central intensity. 

These contours are distorted when the image is taken through a linear polarizing filter. Let 
us use a cartesian coordinate system (p x ,P y ) in the plane of the sky, and assume the polarizing 
filter is oriented along the p y axis. At each point in the image the observed intensity is given as a 
weighted average of the two intensities // and I r , namely, 

I(Px,Py) = h sin 2 9 + I r cos 2 9, (3) 

where 6 is the polar angle of the point relative to the p x axis, so that cos 6 = p x /p and sin 6 = p y /p. 
Using this equation, the contours of the halo were constructed and are shown in the right panel of 
Figure 2, again using a separation of one-half magnitude between the contours. The distortions 
are very evident, even for the innermost contours. 

The detectability of Lya halos around high-redshift galaxies was discussed in detail in Paper 
I. We have found that although difficult, detection of halos at z s ^ 10 might be feasible from 
space. For example, given the upper limit on the cosmic infrared background derived by COBE at 
1.25 /jm (Hauser et al. 1998), one could achieve a signal-to-noise ratio S/N = 10 after 10 hours 
of integration on an 8-meter space telescope (such as the Next Generation Space Telescope) for 
sources at z s ~ 10 which possess a Lya luminosity higher by an order of magnitude than the 
galaxy discovered by Dey et al. (1998). 

Such sources emit N a = 6 x 10 54 s -1 , and might be found in wide-field surveys, based on 
the broad number-flux distribution which is predicted for high-redshift galaxies (see Figure 2 in 
Haiman & Loeb 1998b). 

4. Conclusions 

As shown in paper I, the special type of Lya halos associated with sources surrounded by 
the neutral IGM in a Hubble flow before reionization, can be clearly characterized by their light 
profile and spectrum, which have a universal character. In this Letter we have shown additionally 
that the polarizations associated with these Lya halos are quite large, of order tens of percent, 
and have a particular universal behavior as a function of impact parameter. 

The polarization properties illustrated in Figure 1 could potentially be of importance in 
providing a critical test for distinguishing between Lya halos of the type considered here as 
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opposed to halos due to some other cause. In addition, the polarization signature can improve the 
signal-to- noise ratio in separating faint Lya halos from an unpolarized background light. 

This work was supported in part by the NASA grants NAG5-7768 and NAG5-7039 (for AL). 
The authors gratefully acknowledge helpful conversations with Alex Dalgarno and Kenneth Wood. 
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Fig. 1. — Upper panel: The solid lines give the frequency integrated intensities /, //, and I r versus 
the impact parameter, p. The dashed line is the I found in paper I using isotropic, unpolarized 
scattering. Lower panel: The degree of polarization II versus impact parameter. 
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Fig. 2. — Predicted images of the Lya halo in the plane of the sky, as seen in total light (panel 
a) and as seen through a linear polarizing filter oriented vertically (panel b). The contours are 
separated by one-half magnitude, starting with the central intensity. 



